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The gauge-fixing parameter dependence of two-point gauge variant correla- 
tion functions is studied for QED and QCD. We show that, in three Euclidean 
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fSj ■ dimensions, or for four-dimensional thermal gauge theories, the usual proce- 

dure of getting a general covariant gauge-fixing term by averaging over a class 
of covariant gauge-fixing conditions leads to a nontrivial gauge-fixing parame- 
ter dependence in gauge variant two-point correlation functions (e.g. fermion 



Oh! 

^ • propagators). This nontrivial gauge- fixing parameter dependence modifies the 



large distance behavior of the two-point correlation functions by introducing 
additional exponentially decaying factors. These factors are the origin of the 
gauge dependence encountered in some perturbative evaluations of the damp- 
ing rates and the static chromoelectric screening length in a general covariant 
gauge. To avoid this modification of the long distance behavior introduced by 
performing the average over a class of covariant gauge-fixing conditions, one 
can either choose a vanishing gauge-fixing parameter or apply an unphysical 
infrared cutoff. 
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I. INTRODUCTION 



Physical observables such as thermal damping rates and the Debye screening length 
which are determined by the position of poles (or generally, singularities) in correlation 
functions are gauge invariant quantities However, for a general covariant gauge, a 

gauge-fixing parameter ^ dependence in damping rate was reported |Q in a perturbative 
evaluation of the fermion damping rate. Similar gauge-fixing parameter dependence was 
also encountered in the perturbative calculation of the Debye screening length at the next- 
to- leading order [^. A way to extract the gauge independent damping rate and screening 
length is to introduce an unphysical infrared cutoff f^. In this article, we examine the 
gauge-fixing parameter ^ dependence of two-point correlation functions. As is well known, 
a conventional way of getting a general covariant gauge-fixing term {dAy'/{2S,) involves 
performing an average over a class of covariant gauge conditions. We shall show that, in 
three Euclidean dimensions or for four-dimensional thermal gauge field theories, this average 
over different gauge conditions generates a nontrivial ^ dependence in two-point correlation 
functions of gauge variant operators such as the fermion propagator. Specifically, for QED 
and to leading order of QCD, this gauge-fixing parameter dependence alters the long range 
behavior of the two-point correlation functions by an extra exponentially decaying factor 
with the exponent depending on ^. This is the origin of the ^ dependence encountered in 
perturbative calculations of the damping rate and the Debye screening length This 
gauge-fixing parameter ^ dependence is an artificial fact due to doing the average over a class 
of gauge conditions since this average includes gauge conditions dA = f with / containing 
long wavelength fluctuations. Choosing the Landau gauge ^ = which means taking the 
"no average" limit removes the modification of the long distance behavior produced by 
averaging over gauge conditions. Another way is to introduce an unphysical infrared cutoff 
as suggested in reference to suppress the contributions to the gauge condition average 
from the infrared fluctuations. If the physical content of the theory is gauge invariant in 
the sense that the gauge constrain dA = is equivalent to other physical gauge constrains 
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when quantizing the theory, choosing ^ = then yields the correct physics ( e. g. correlation 
length) . 

In next section, we shall first study generally how the average over a class of gauge choices 
produces the ^ dependence of the propagators of charged particles in QED where it is com- 
pletely solvable. The ^ dependence of the propagator of a charged particle may be expressed 
simply as a ^-dependent multiplicative factor which becomes an exponentially decaying fac- 
tor for large spacetime argument in three Euclidean dimensions or four-dimensional thermal 
QED. Thus, if one extracts the correlation length or the damping rate from the propagator 
of a charged particle, a ^-dependent correlation length or damping rate is obtained. We 
then explain why an unphysical infrared cutoff can get rid of this ^ dependence and produce 
the same result as the choice ^ = 0. In section III, we perform parallel analysis for QCD. 
Since a complete solution could not be achieved, we only do a leading order perturbative 
calculation. To the leading order, the gauge-fixing parameter dependence for the fermionic 
propagator is the same as in QED with an effective charge. The ^ dependence of the static 
(chromo)electric screening length is discussed. We draw conclusions in section IV. In Ap- 
pendix A, details about the evaluation of a thermal integral is given. In Appendix B, we 
present another way of deriving the ^ dependence in fermion propagators based on Ward 
identities for the proper vertices. 

II. ^ DEPENDENCE OF A CHARGED PARTICLE PROPAGATOR IN QED 

A. A functional derivation 

In this subsection, we shall derive the gauge-fixing parameter ^ dependence of the prop- 
agator of a charged particle in QED in a general covariant gauge. Since the derivation does 
not depend on whether the charged particle is a scalar or a fermion, we focus on the fermion 
case. Throughout this paper, we work in Euclidean spacetime. Real time results can then 
be obtained by analytic continuation. 
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The fermion propagatoiQis defined by the Euchdean functional integral as 

Gf{x,y) = {il}{x)i){y)) 

\VA][V^[V,P] exp (- / £e j ^{x)ij{y)5{dA - f) , 



where we have introduced the gauge-fixing condition 

dA{x) = fix) 

with f{x) being an arbitrary function. Denoting the charge of the fermion by q, it 
hard to justify that the gauge transformation change of variables 

^{x) ^'{x) = e''''^^''^{x) 
A{x) A'{x) = A{x) - dX{x) , 

with A satisfying 

d'X{x) = fix) , 

changes the gauge-fixing condition 

dA = f dA' = 

and thus gives 

Gf{x,y) = Gf=o{x,y)e''^'^^^-'^y^K 
Defining Green's function A{x — x') by 

5^ A(x — x') = S{x — x) 

enables us to write 



^We use (• • •) to represent the Euclidean time ordered product. 



X{x) = J dzA{x-z)f{z). 



(2.8) 



The usual procedure of getting the general covariant gauge now involves averaging over f{x) 
with a weighting factor exp{ — ^ / dxp}. Let us denote the fermion propagator in a general 
covariant gauge by G^{x,y). Then 



G^{x,y) = J [I?/]GKx,y)exp|-l J dzf{z) 
= |[P/]exp|- J dz 
where we have defined the linear source J{z; x, y) as 



^P{z)+iJ{z-x,y)f{z) 



Gf=o{x,y) 



(2.9) 



J{z] X, y) = — g [A{z — x) — A{z 



A straight forward evaluation of the gaussian functional integral in Eq. (|2.9[) yields 



(2.10) 



G^{x,y) = exp i — — J dz [A( 



z-x) - A{z - y)f \ Gf=o{x, y) 



where the normalization factor for the functional integral has been chosen so that 



(2.11) 



[I?/]exp 



2^ 



dzf{z) 



(2.12) 



It is trivial to check that the Landau gauge choice ^ = corresponds to the gauge choice 
/ = 0. This is expected since the weighting functional exp{ — ^ / dxf"^} for the average 
allows / to fluctuate around / = with the variance proportional to ^. Setting ^ = 
confines / to be 0. The gauge-fixing parameter dependence appears as a multiplicative 
factor. Result ( p.ll|) is also valid for the propagator of a scalar charged particle with charge 
q. Equation (|2.11| ) has been derived by other methods and discussed for four- dimensional 
QED §0. 

By including the photon source and more pairs of the Green's function A{x) in the 
source ( |2.10| ), Eq. ( |2.11| ) may be generalized to cases where the correlation functions contain 
external photon lines and additional fermion lines. The / functional integral is still a gaussian 
integral. We omit the algebraically complicated intermediate steps which are completely 



parallel to those for the fermion propagator. The ^ dependence for a correlation function 
with 2n external fermion legs is 

= exp \^-^ j dz X: {^{z-x,) - ^{z-y,)) + j dz'A{z-z') d,,-j{z') 

X (^(xi)V^(x2) ■ ■ -^(xO V^(yi)V^(i/2) ■ ■ ■ ^(j/n) e^^/^'-^)/=o , (2.13) 

where we have introduced the photon source term iq J j ■ A. Taking derivatives with respect 
to the photon source j gives insertions of the photon fields in the correlation function. We 
note that the ^ dependence is totally factorized. Therefore, the choice ^ = is equivalent 
to the gauge choice / = 0. 

B. ^ dependence in four and three dimensional QED 

We now study the behavior of the multiplicative ^-dependent factor in four and three 
dimension spacetime. We do not consider the case d = 2 since there infrared divergences 
are so serious that the charged particles are confined 0. To facilitate the notation, let us 
define 

I{x) = ^Jdz[A{x-z)-A{z)f (2.14) 

so that 

G^ix,y) = e-^^''^^-y^Gj=oix,y) . (2.15) 



In momentum space, I{x) can be expressed^ as 



I[x) = - f {dk) -^-(2 - e^'^^ - e-'^"") = f {dk)^{l - e*'^^) , (2.16) 



^Here, we have made a change of variable k — > —k. In three dimensions, this causes an infrared 
problem which does not really matter if we only use regulators invariant under the inversion of 
momentum. 



where / (dk) represents the appropriate momentum integral conjugate to the spacetime. In 
d dimensional Euclidean spacetime, 



I{x) 



d'^k 1 



,1 - e 



ikx\ 



(2.17) 



which can be evaluated as 



r ds s'-"/' exp 
{iny/^ Jo 

r(d/2 -2) ^ 2\ 2-^/2 



X 



{Any 



Setting d = 4, Eq. ( |2.15| ) takes the form 



Gd^,y) = exp 



\(47r)2d-4 



(2.18) 
(2.19) 



where we have introduced the renormalization scale /i by 

q'^ = qURf^'-' (2.20) 

with q^^ and being the renormalized charge and gauge-fixing parameter. 7 in Eq. (p.l9|) 
is the Euler constant. The ultraviolet divergence appearing in the gauge dependent factor 
gives the usual gauge dependence of the fermion wave function renormalization factor in 
agreement with previous results P,PJTI1|] as is the (x — y)^ power modification factor [Q. 
In four dimensions, the ^ dependence in the fermion propagator modifies the long distance 
power law behavior. Due to massless photons, the fermion propagator does not exhibit a 
pole in momentum space but a branch cut with the behavior (p^ -|- m2)^(^+'^) [0!0]- The 
^-dependent modification factor (|2.19|) leads to a ^-dependent z/ and thus a ^-dependent 
on-shell condition for fermions in four dimensional QED 0. 

We now consider QED in d=3 Euclidean spaceQ. We can view this as a dimensionally 
reduced field theory of a four dimensional scalar QED at the high temperature limit |[19|,p0|. 



^For a real three dimensional QED, there is evidence [|l2 14] showing that the charged particles may 
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I{x) is both ultraviolet and infrared finite. Eq. (|2.18|) reads explicitly /(x) = |x|/(87r). We 
have then 

G^(x,y) = e-^l^-ylG;=o(x,y). (2.21) 

Hence the usual average over covariant gauge conditions has introduced a ^-dependent ex- 
ponentially decaying factor to the propagators of charged particles. Thus, the correlation 
length also acquires ^ dependence. If we extract the correlation length from a gauge variant 
two-point correlation function, different values give different answers. Choosing ^ = gives 
the answer corresponding to the original gauge theory quantized by the gauge condition^ 
dA = 0. This justifies a previous claim on the preference of the Landau gauge choice ||15|| . 
We note here that if only the gauge invariant correlation functions are considered, there are 
no ^ dependences inside the correlation functions which states that all ^'s are equivalent for 
gauge invariant correlation functions. 

We now examine the origin of the decaying factor in Eq. (|2.21|) and explain why an 



unphysical infrared cutoff can remove this artificial ^-dependence. In previous subsection, 
we found the relation 



Gf{x,y) = exp ^-tq J dz [A{x - z) - A{y - z)] f{z)jGf=o{x,y) (2.22) 
which basically states that different choices of / are equivalent. However, this statement is 



be confined in three dimension QED due to infrared divergences. Of course, one can consistently 
add to the theory a topological mass term for the photon field without breaking the local gauge 
symmetry [1^.1^. This topological mass term can be even dynamically generated by interacting 



with the fermions [15-18|. With this topological mass term, the charged particles are no longer 



confined |14,15|. 

^It should be mentioned that if we view the three dimension theory as the high temperature limit 
of a four dimension theory, the gauge condition dA = really corresponds to the Coulomb gauge 
choice in the original four dimensional theory. 
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based on the assumption that the factor 

exp J dz [A{x - z) - A{y - z)] f{z)^ (2.23) 

does not alter the long distance behavior of the propagator or, in particular, the position 
of the physical pole appearing in the Fourier transform of Gf{x,y). If f{z) is a localized 
function so that the exponent in factor ( |2.23| ) vanishes as x or y becomes large, Gf{x,y) 
has the same large distance behavior as Gf=o{x,y). We can then conclude that the gauge 
choice dA = / is equivalent to the choice dA = 0. However, when performing the average of 
gauge conditions, nonlocalized /'s are not excluded. It is not hard to see that the weighting 
functional^ 

exp{-l/d.f(.)) =exp{-i/|l /»/(-*)) (2.24) 

contains long wavelength modes {f{k) with wavelength 1/k longer than the separation be- 
tween X and y). It is the inclusion of these long wave modes in the average that yields an 
exponentially decaying factor depending on |x— y|. For the gauge conditions with / having 
wavelength shorter than |x— y|, the x,y dependence in phase factor ( p.23| ) is washed out 
after summing over different short wavelength contributions. In another word, the short 
wave fluctuations do not suffice to change the behavior of the long range correlation after 
the average. Indeed, I{x) does get its main contribution from the infrared region with k 
being order l/|x| or less. This generates a piece proportional to |x| and therefore results 
an exponentially decaying factor. Employing an unphysical infrared cutoff to suppress the 
contributions from these long wave /'s can eliminate the gauge-fixing parameter dependence 
in the correlation length. Explicitly, introducing an infrared cutoff fcmin to integral /(x), we 
have, as |x| oo. 



fc>fcmin (27r)3 k"^ 27r2 Jk^^^ k"^ 



sin 


A; X 


k 


X 



oo 



sm s 



x| . , , ds—r - (2.25) 



-"Here f{k) is the Fourier transform of f{z). 



9 



which does not depend on x. Therefore, the average over the short wave /'s does not 
modify the long distance behavior of the propagators but an overall constant. This explains 
why an unphysical infrared cutoff proposed in reference can remove the ^-dependent 
modification of the long distance behavior. 

Since lower dimension field theories are more sensitive to the infrared region, this expo- 
nentially decaying factor does not appear ind = 4. In four dimensions, only the power law of 
the propagator is changed. For four-dimensional thermal field theories, the imaginary time 
formalism leads to a Euclidean functional representation for thermal correlation functions 
with one dimension of the spacetime compactified. Therefore, we expect that the average 
over different covariant gauge choices may also cause serious modifications to the two-point 
correlation functions. 



C. Four-dimensional thermal QED 

We now study the four- dimensional thermal QED which is equivalent to QED in three 
spatial dimensions plus an additional compactified imaginary time dimension. As such, 
the momentum in this dimension is discretized so that the momentum integral for I{x) in 
Eq. ( p.l6|) is a "sum-integral": 

ko 



(2.26) 



Since there is the zero temperature part contributing to the sum-integral, it is ultraviolet 
divergent. Using the dimensional regularization, we calculate this sum-integral in Appendix 
A and simply quote the result ( |A1|) as: 



/(r,x) 



(47r) 



r(4vrT 



2\d/2~2 



A-d 
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+ 



(47r 



^ In f 1 + e-^"^l"l - 2 cos(27rTr)e-2-^l"l) + ^ + 0(ci - 4) . (2.27) 
n)'^ ^ ^ Stt 



Here I{x) has been written as /(r, x) with r and x being the imaginary time and spatial 
coordinate respectively. It is not hard to check that for large |x| or high T 
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/(r,x)^^ (2.28) 

which is in agreement with the analysis for the three dimensional Euclidean theory dis- 
cussed in previous section. Therefore, the correlation length contains a ^-dependent piece 
eg2T|x|/(87r). 

So far the correlation functions studied are all defined in imaginary time. To study the 
damping rate, a retarded real-time correlation function is required. We can analytically 
continuate it into real-time by replacing r with it to obtain the corresponding real-time 
propagator .0 Performing this analytic continuation for /(r, x) gives 

2 



/(r,x) = ^(4.T^r/-^ 



4-d 
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— In f 1 + e~^"^l"l - 2 cosh(27rTt)e-2'^^l"l) + + 0{d - 4) . (2.29) 



(47r)2 V ' ' ^ Stt 

For large |x| or we have 

J(t,x) ^^[\t\ e{\t\ - |x|) + |x|^(|x| - m . (2.30) 

Inserting this large coordinate argument behavior into Eq. ( |2.15D gives 



G'^(t,x) ~ exp 1-^^ [\t\ 9i\t\ - |x|) + |x|^^(|x| - |t|)]| G^=o(t, x) , (2.31) 



in 



where we have shortened]]] the Green's function notation G{x, y) to G{x—y) because of the 
spacetime translation invariance. The ^-dependent large time damping factor causes a ^- 
dependent damping rate if we use the general covariant gauge. To avoid any modification of 
the large time behavior coming from averaging the gauge conditions, we can choose ^ = 0. 



"This analytic continuation yields real-time ordered correlation functions. However, it is not hard 
to show that the exponentially decaying factor we are concerned is the same as that for the retarded 
correlation function. 

^We have implicitly switched back and forth between the notations x and (t, x) for the spacetime 
coordinate. 
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This choice yields the damping rate for the original gauge theory, if the theory is invariant 
for different gauge constrains (without involving any average) when quantizing it. We also 
provide a more familiar derivation of Eq. based on the proper vertex Ward identities 



in Appendix B where Equation ( [B5| ) shows that the ^ dependence we discussed above is the 
same ^ dependence as reported in reference 0]. We now examine again why an unphysical 
infrared cutoff can also remove this ^ dependence f^. Introducing a small mass term to 
cutoff the k integral for /(r, x) at the infrared region, and using the usual contour trick to 
do the ko sum, the finite temperature part of I{t,x) is expressed as 



cim2 [J (27r)3 2V^Tm^ 
where ni^u) is the Bose distribution factor 



^(^) - ^^3^ (2-33) 

with j3 = 1/T. Doing the analytic continuation r it, we find that as t or |x| becomes large, 
the Riemann-Lebesgue lemma kills t and |x|-dependent part in I^'^\t,x). Hence the large 
coordinate argument behavior does not obtain any ^ dependence but an overall constant. 
This shows why an unphysical cutoff in reference can remove the gauge-fixing parameter 
dependence in the damping rate. 



III. GAUGE-FIXING PARAMETER DEPENDENCE IN THERMAL QCD 

The gauge-fixing parameter dependence we studied for QED can also be derived by using 
Ward identities. Since to leading order in QCD, Ward identities are the same as that in 
QED, we expect that at the leading order, we should be able to find similar gauge-fixing 
parameter dependence as occurs in QED. This is indeed the case. Since the derivation for 
QED is completely parallel to the leading order derivation for QCD, it suffices to just provide 
the derivation for QCD at the leading order. 
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A. fermionic damping rate 



Let us consider the fermion propagators. To simplify the notation, we shall suppress 
the color indices of the fermion fields. We shall use the Ward identities to derive the result 
for the QCD fermion propagator analogous to the result ( ^.11 ). Taking the derivative with 
respect to ^ which introduces an insertion of the gauge-fixing term in the functional integral, 
we obtain 



d 



1 



dz{^p{x)^{y)[dA^{z)f) 



(3.1) 



where ip a^nd A represent the fermion fields and the gauge fields respectively and a is the 
color index in the adjoint representation. Using the equation of motion for the ghost field, 
the BRS transform of the correlation function {il!{x)ilj{y)dA°'{z)c"'{z)) produces the relation 



{^P{x)ij{y)[dA-{z)f) = tg^{tP{x) c\x) - c\y) T'^Piy) [dA^{z)]c^{z)) , 



(3.2) 



where T are the generators in the fermion representation and and c^{z) are the 

Faddeev-Popov ghost fields. Noting that 



{c\x)A'^iy)c'^iz))=0. 



(3.3) 



we have, at the leading order. 



{^{xMy) [dA\z)f) =-zgCmx)T^{y) [dA^z)]) [{c\x)(^{z)) - {c\y)c^{z))] + 0{g') 

= -tg^{ijix)T^i;iy) [dA'^iz)]) [A,^ix-z)-A,^{y-z)] + 0{g') , (3.4) 



where the ghost propagator Agh(a;) is defined as 



{c\x)c\y)) ^ 5^' A,^{x ~ y) 



Similarly, the BRS transform of {il){x)T""il){y)c°'{z)) yields 



(3.5) 



(^(x)TXi/) [dA\z)]) = -igi{i,{x)i:''T''i^{y)) W^{x-z) - \^{y-z)] + 0{g') 

= -tg^CFG^{x,y) [A,^{x-z) - A,^{y~z)] + 0{g^) , (3.6) 



13 



where Cp is the Casimir for the fermion representation. Combining results (|3.4| ) and ( p.6|) 
above gives 

^^G^{x, y) -^/Cf / dz [\^{x-z) - Agh(2/-^)]' G^{x, y) + 0{g') . (3.7) 
Perturbatively, 

Agh(a:) = A(x) + 0(/). (3.8) 

Therefore, to the leading order, the above differential equation shows that the leading ^ 
dependence of the fermion propagator can be expressed as 

G^{x,y) ~exp|-^^^/(x-?/)|G/=o(a;,l/). (3.9) 

We note here that a QED derivation may be obtained simply by ignoring all the color 
indices, changing T'' — * 1, (7 — g, and setting all the leading approximation to be exact. 
For QCD, the leading gauge-fixing parameter dependence of the fermion propagator is 
the same as the exact ^ dependence of a fermion propagator in QED with effective charge 

1 /2 

gG-p . As discussed before, at finite temperature, if we extract the damping rate of fermionic 
excitations from the fermion propagator calculated in a general covariant gauge, we shall 
get a ^-dependent damping rate. This is the dependence reported in reference 0. We can 
either choose Landau gauge ^ = or use an unphysical infrared cutoff to get rid of the 
modification of the long time behavior due to taking the average over gauge conditions. 

B. Static electric screening length 

We now turn to consider the static chromoelectric screening length to next-to-leading 
order. Since the relevant energy scale is gT, it is convenient to use the dimensionally 



This is also true after performing the Braaten-Pisarski resummation since the resumed proper 
vertices and propagators still satisfy the QED-type Ward identities . 
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reduced effective tlieory (so called EQCD) which involves only the static gauge fields to 
describe thermal QCD |2^. To the order we are concerned, we only need to study a three 
dimensional Euclidean effective theory with Lagrangian 

C^QCB = \f^F^ + \{D,A,r{D,A,r + (3.10) 

where the covariant derivative Di is defined as 

{D,Aor = d.Al-igr'^A\Al (3.11) 

and mei is the leading order Debye mass with the value 

ml = \{c^ + ^)g'T\ (3.12) 

Here, Ca is the Casimir of the adjoint representation for the gauge group and A^f is the 
number of fermion flavors. Of course, we shall study, for a general covariant gauge, the 
gauge-fixing parameter ^ dependence of the static propagator D"''(x, y) of the Aq fields. We 
have not explicitly included the gauge-fixing term nor the ghost fields term. We employ 
the same steps as in the last subsection for calculating the gauge dependence of a fermion 
propagator. All we need to do is to change all the fermion fields into the Aq fields and 
replace the generators T" by the generators in the adjoint representation under which Aq 
transforms. Thus, the effective charge squared g'^C-p appearing in Eq. (|3.9|) is changed to 
g^Cp^. After these replacements, we obtain 

(x,y) ^ exp |-i^|x- y|| Df^o(x,y) . (3.13) 

Performing the Fourier transform gives again a ^-dependent singularity in the propagator. 
This dependence appears in the result^ found in references |^J^ where the ^ dependence 

^There, a ^-dependent term ^g^A^mciT/(47r) contributes to the self energy which indicates a cor- 
rection ^g^NT / {Stt) to the leading screening mass rriei. N is the number of colors, or equivalently, 
the Casimir of the adjoint representation for SU(N) gauge group. 
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is removed by introducing an unphysical cutoff at the infrared region of the loop integral 
involved in the self energy evaluation. 

To avoid confusion, we like to add following comment. It appears that in reference [Q, 
the Feynman gauge choice, ^ = 1, coincides with the result obtained by introducing an 
unphysical infrared cutoff. The free gauge boson propagator 

= + (3.14) 

contains the combination ^ — 1 as the coefficient of PiPj/p'^- At the one- loop order, it 
can be shown explicitly that the part in the gauge boson propagator proportional to 
PiPj/p'^ contributes to the self energy and causes a shift of the position of the pole of Aq 
propagator. Introducing a technical infrared cutoff removes this contribution to the pole 
position of this propagator from the PiPj/p^ part of the gauge boson propagator. This 
can be mimicked by choosing ,^ = 1. However, the transverse part of the gauge boson 
propagator receives additional radiative corrections while the pure longitudinal part does 
not as a consequence of the Ward identity. Due to these radiative corrections, this transverse 
piece becomes less singular at the infrared region PJ20[] and a physical infrared cutoff gets 
induced. Consequently, its piPj piece does not shift the position of the pole in Aq propagator. 
On the other hand, the longitudinal part remains the same and does shift the position of 
pole except in Landau gauge with ^ = 0. Thus, the Landau gauge choice produces the same 
result as a naive choice of Feynman gauge where in addition the radiative corrections to the 
gauge boson propagator are ignored. For theories containing a topological mass term, the 
gauge boson propagator has the form 

p^ + "^topo [p^ + "^topo) p^ r 

where mtopo is the topological mass. Here it is easy to see that the infrared behavior of 
the longitudinal part is different from that of the transverse part. Choosing a vanishing 
gauge-fixing parameter is the same as putting an infrared cutoff. 
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IV. CONCLUSIONS 



In conclusion, we have shown that in three Euchdean dimensions or for four-dimensional 
thermal gauge field theories, the usual averaging procedure for getting a general covariant 
gauge-fixing term may introduce gauge-fixing parameter dependent modifications to the 
large distance behavior of the gauge dependent correlation functions. The gauge-fixing 
parameter dependent modification to the large distance behavior of the correlation function 
is the origin of the gauge dependence encountered in some perturbative evaluations of the 
damping rate and the Debye screening length. Choosing a vanishing gauge-fixing parameter 
(Landau gauge) or introducing an unphysical infrared cutoff enables us to avoid this gauge- 
fixing parameter dependent modification at the long distance introduced by averaging over 
a class of gauge conditions. If the theory is gauge invariant in the way that it can be 
quantized by using different gauge constrains (without involving any average), we can then 
extract physics from gauge variant propagators evaluated in a general covariant gauge by 
choosing a vanishing gauge-fixing parameter. 
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APPENDIX A: EVALUATION OF I{X) FOR THERMAL GAUGE THEORY 



For a compactified spacetime, the spacetime point x is understood as a; = (r, x) . Using 
the dimensional regularization to regulate the ultraviolet divergence, we evaluate the integral 
I{x) defined by Eq. ( |2.1(j| ) as 

d'^-^k 1 



1 - e 

(27r)'^-i (k2 + A;2)2 ^ 

(2vr)^-i (k2 + klY 



(27r)3 (k2 + klf 



+T 



d'^-^k 1 
(27r)^-iF' 



l-e*'^-'') + 0(d-4) 



Tlx! 



(47r)(rf-i)/2 



57r 



0(rf-4) 



(4 



4-d 



+ 7 



(An) 



(47r)' 



ln(l - e-2'^^l''l+*2'rTr-j _^ _ g-27rT|x|-i27rTr^ 

2 



+ ZM + o(rf-4) 



57r 



4-rf 



+ 7 



(47r) 



^ In f 1 + e-^-^l^l - 2 cos(27rrr)e-2-^l^l) + ^ + 0(rf - 4) . (Al) 
7r)2 ^ ^ Svr 



APPENDIX B: ALTERNATE DERIVATION BASED ON THE PROPER 

VERTICES WARD IDENTITIES 

Consider the fermion self energy diagrams in QED. The only place where the ^ parameter 
can enter is in the longitudinal part of the photon propagator E^k^k^j/k'^. Therefore, the 
derivative of the self energy with respect to ^ is 

^S(p) = g2 I {dk) '^T,{k, p,p + k)G{p + k)T,{-k,p + k,p) 

+^l{dk)^T,^{k,-k,p,p), (Bl) 

where T ^{k^p^p + k) and T^y{k, k' ,p,p + k + k') are the proper photon-fermion three-point 
and four-point vertices respectively. For the proper vertices, our convention is that the last 
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two momentum arguments are the momentum of fermion legs while the k and k' denote the 
photon momenta. Ward identities give 



k^T,{k,p,p + k) = G-\p + k)- G-\p) 
k^k,T^,{k, -k,p,p) = G-\p + k) + G-\p -k)- 2G-\p) . 



Inserting these identities into Eq. (Bl) producesj^ 



^S(p) = -q'G'\p) I (dk)^ [G{p + k)G-\p) - 1 



Since G{p)G ^{p) = 1, taking the derivative with respective to ^ yields 

^Sb) = ^^G-\p) = -G-\p)j-^G{p)G-\p) . 
Combining the two equations above, we obtain 

j-^G{p) = q' j{dk)^^[G{p + k)-G{p)] , 
which can be written in coordinate space as 

^G(x) = g^/(rfA:)l(e^'^^-l)G(x). 
It is straight forward to solve the differential equation above to get 
G^{x) = exp j idk)^il - e^^^^)} G5=o(x) . 

This is the result derived by functional methods in the main text. 



(B2) 



(B3) 



(B4) 



(B5) 



(B6) 



(B7) 



^'^This equation agrees with the equations appearing in references [^J5| where the gauge-fixing 
parameter dependence was examined. 
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